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INTRODUCTION 


The tasks proposed in our last proposal have been met and additional 
information generated. The July 1986 Progress Report has been modified and 
sent to The Journal of Crystal Growth to be published. The article, entitled 
"The Effect of Gravitational Field Strength on Steady Thermal Convection in 
Bridgman Stockbarger Crystal Growth," is an extension of the data presented at 
The Gordon Conference in 1985. These results should be of considerable 
interest to the crystal growth community. 


DISCUSSION OF PROGRESS 


The primary focus during this reporting period has been on completing a 
simulation of thermos olutal convection. The configuration employed was based 
on seme MIT experiments using a Germanium-Silicon charge in a quartz ampoule 
with an aspect ratio of 8 and a rather large adiabatic zone as given in Fig. 
1. Two thermally stable cases were investigated: case A, solutally unstable, 
and case B, solutally stable. Both use the same equilibrium phase 
diagram(Fig. 2) and thermophysical properties except the sign of the solutal 
coefficient of volume expansion is reversed. Initially the ampoule is placed 
midway in the adiabatic zone and is in a steady state thermal convection mode 


with a zero pull rate(V 
interface temperature(T f = 


p = o) - 
937. A°c) 
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The gravity level is 9.8(10 ) m/ s and an 

corresponding to a zero concentration in 


the melt presides. 


At time equal to zero the ampoule pull rate is increased to 1.39(10 ) 

m/s and a uniform melt mass fraction, C, of 0.2 is impressed on the system. 
Case A will be described first followed by differences between the cases. The 
interface position as the solidification progresses through the ampoule is 
shown in Fig. 3. No movement is indicated in the first 200s as the interface 
temperature adjusts to its new lower equilibrium value. The crystal was not 
allowed to melt. A more realistic situation would have been a melt back 
followed by a resolidification. Between 200s and 700s a transient growth 
phase where interface growth rate, V^, and pull rate are not the same. The 

initially planar interface now becomes concave and the curvature continues to 
change with time. Following this a steady state phase from 700s to 1370s 
occurs in which pull rate and growth rate are equal. Bulk concentration 
continually changes during all phases of the process. Finally, the side of 
the ampoule is lowered to a position within the adiabatic zone, and since the 
top is insulated, rapid solidification ensues until the entire charge is 
crystalized at 1810s. 

Figure 4 shows the concentration field throughout the charge at 748s. 
Notice that the interface is not at a constant axial location but varies with 

radius. In this plot it lies between z = 4.05 and 4.15. Inspection will 

reveal that the interface is concave. The liquid region immediately in front 

of the interface indicates a thin layer (approximately 0.5 cm) of rapidly 

decreasing concentration. This is the concentration "boundary layer". There 
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is also a radial gradient in the liquid(say at z = 5.5) outside this layer due 
to advection(a single cell, more later) in the bulk. This contrasts with the 
dual celled solutally stable case B which has a corresponding smaller 
gradient . 

At 1810s the charge is completely solidified. The C distribution in the 
solid will be examined using Fig. 5 and Fig. 3. 

3.5 < z < 3.75, (0 < t < 450s) 

Initial transient period, interface adjusting to the step change in bulk 

concentration. The interface temperature must drop from 937 to 918 °c. 

= 0 until about t = 200s. Then increases from zero to sane value 

greater that and Cs increases. Cs increases because solute cannot 

"escape" to the bulk but instead is "trapped" in the solid. 

3.75 < z < 4.0, (450 < t < 700s) 

V^_ decreases and then increases. Corresponding to this Cs increases and 

then decreases. This always happens and others have noticed it [e.g., 
Snith, Tiller, Rutter, Canadian J. of Physics v33 no. 1 (1955) 723-745, 

and Bourret, Derby, Brown; J. Crystal Growth 71 (1985) 587-596.]. 

4.0 < z < 4.9, (700 < t < 1400s) 

"Steady state" growth with continuous change in 

1.39(10 ^) m Is = slope of curve in Fig. 3. Note 
segregation, 

4.9 < z < 7.5, (1400 < t < 1800s) 

Side of the ampoule now entirely within the adiabatic zone. No energy 
enters the charge since the top of the ampoule is insulated. exceeds 

V as shown. This increase in dictates an increase in C as previously 

mentioned, but the overall species concentration demands a general 
decrease. The characteristic increase in C in the vicinity of the "last 
to freeze" is not quite as dramatic as it should be because of the rather 
large computational volumes and the averaging used. 

Two curves corresponding to those just presented are given in Figs. 6 and 
7 for solutally stable case B. Differences between the two cases will be 
noted below with the aid of Table 1. Case A is primarily single celled 
advection with a second, much weaker cell appearing for short time intervals 
but always fading. The melt velocity increases from the initial thermally 

-2) 

driven condition to 3.62(10 cm/s and then decreases. The interface is 
essentially planar from the ampoule centerline to about r/R =0.5 after which 
it becomes concave. 

In case B, two counterrotating cells are present for most of the process 
except at the very beginning and end. The melt velocity is maximum at 
-2 

2.0(10 ) cm/s at time equal to zero and then steadly decreases to zero. The 


C-bulk. V r = V = 
f p 

the large radial 
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upper cell is always stronger than the lower cell but is extinguished at the 
end of the process before the lower cell. The magnitude of the maximum 
velocity is always less than case A at corresponding times. The concave 
interface is usually planar between r/R = 0.0 and 0.7. 
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Solutally Unstable 





Specific heat 
0.39 Cal/ g K 
Prandtl number 
5.88 CIO -3 ) 

Radius R 

0.9525 cm 
T c = 765°C 
Tf = 937°C 
T h = 107 0°C 

Thermal coefficient of volume expansion 
5 (10' 4 )/K 

Melt thermal diffusivity 
0.187 cm 2 / s 

Crystal thermal diffusivity 
0.0815 cm 2 / s 

Wall thermal diffusivity 
0.1247 cm 2 / s 
Density 

5.48 g/cm 3 
Kinematic viscosity 
1.10 (10‘ 3) cm 2 / s 
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Figure 1. Schematic of Bridgman-S tockbarger System and Thermophysical Properties 



Figure 2. Phase Diagram 



SYMBOL 



CO 


convec t ion . 
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z along the ampoule Can) 

Figure 4. Charge concentration field when partially solidified at 506 seconds. 
Thermally stable, solutally unstable. 
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Figure 5. Charge concentration field when completely solidified at 1810 seconds. 
Thermally stable, solutally unstable. 
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Figure 6. Interface position versus time for thermally stable, solutally 
stable convection . 
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z along the ampoule Can) 

Figure 7. Charge concentration field when completely solidified at 1800 seconds. 
Thermally stable, solutally stable. 



